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ORIGINAL ARTICLE 

Quantitative Systems Pharmacology Model of NO 
Metabolome and Methemoglobin Following Long-Term 
Infusion of Sodium Nitrite in Humans 

K Vega-Villa 1 , R Pluta 2 , R Lonser 23 and S Woo 1 

A long-term sodium nitrite infusion is intended for the treatment of vascular disorders. Phase I data demonstrated a significant 
nonlinear dose-exposure-toxicity relationship within the therapeutic dosage range. This study aims to develop a quantitative 
systems pharmacology model characterizing nitric oxide (NO) metabolome (NO A/0 2 " o A/0 3 ") and methemoglobin after 
sodium nitrite infusion. Nitrite, nitrate, and methemoglobin concentration-time profiles in plasma and RBC were used for model 
development. Following intravenous sodium nitrite administration, nitrite undergoes conversion in RBC and tissue. Nitrite 
sequestered by RBC interacts more extensively with deoxyhemoglobin, which contributes greatly to methemoglobin formation. 
Methemoglobin is formed less-than-proportionally at higher nitrite doses as characterized with facilitated methemoglobin 
removal. Nitrate-to-nitrite reduction occurs in tissue and via entero-salivary recirculation. The less-than-proportional increase 
in nitrite and nitrate exposure at higher nitrite doses is modeled with a dose-dependent increase in clearance. The model 
provides direct insight into NO metabolome disposition and is valuable for nitrite dosing selection in clinical trials. 
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The effect of nitric oxide (NO) as a vasodilator agent has 
been studied in vascular homeostasis and pathology. NO is 
synthesized in the arterial vascular endothelium, undergoes 
oxidative and reductive reactions, and yields endogenous 
metabolites: nitrite (NO~) and nitrate (/VO~). 1 Nitrite and 
nitrate serve as NO reservoirs in blood and transport it for 
delivery to various tissues and organs 23 due to NO's short 
half-life in blood (< 2 ms) and limited diffusional capability. 4 
Under normal physiological conditions, NO concentrations 
range from 0.01 nM to 1 uM. 4 In contrast, endogenous nitrite 
and nitrate are found in the nanomolar and millimolar range, 5 
which facilitates their role as NO donors. Unlike therapeutic 
nitrates (e.g., nitroglycerin) that may be associated with tol- 
erance, 6 administration of therapeutic nitrite does not show 
tolerance and has been proposed for the treatment of car- 
diovascular diseases. 7 Preclinical and clinical studies have 
shown that nitrite administration produces vasodilation, cyto- 
protection after ischemia reperfusion injury, and blood pres- 
sure reduction. 1 

Two clinical studies have evaluated the safety of nitrite 
as a NO donor for human administration by short-term infu- 
sion. 89 Nitrite readily diffuses into RBC and interacts with 
hemoglobin; 10 thus, the characterization of nitrite uptake 
and release by RBC plays an important role not only in 
describing nitrite disposition, but also in understanding 
the pharmacodynamics of RBC as a vehicle to transport 
nitrite. Hemoglobin-nitrite interactions form NO, nitrate, 
and methemoglobin 11 (Figure 1). Methemoglobin, formed 
by hemoglobin oxidation, does not bind oxygen due to the 
ferric state of iron in the heme molecule. Under normal 
physiological conditions, methemoglobin is present in trace 



amounts (< 1%); however, injury or oxidizing agents includ- 
ing nitrite may convert hemoglobin into methemoglobin, 
which increases methemoglobin concentrations. Therefore, 
a methemoglobin level above 5% is used as a marker for 
nitrite-related toxicity 912 . 

Hon etal. 9 characterized the kinetics of nitrite, nitrate, met- 
hemoglobin, and iron-nitrosyl hemoglobin (HbNO) in plasma 
and RBC in normal volunteers following intra-arterial short- 
term infusion of sodium nitrite incrementing from 42 to 6,600 
ug/kg/h at 5min intervals, and provided valuable insights into 
nitrite-hemoglobin interactions. Nitrite in plasma declined 
biexponentially, and returned back to baseline concentra- 
tions 3h after infusion cessation. 89 Concentrations of nitrate 
declined slower than those of nitrite. The methemoglobin rise 
was lower than expected at higher nitrite doses, which sug- 
gested nonlinearity in methemoglobin dynamics, and was 
described using two rate constants, 9 limiting the model appli- 
cation to studies with alternative dosing regimens. 

Because therapeutic application of sodium nitrite will 
require prolonged infusions, a Phase I trial was performed 
to assess the safety and PK of increasing doses of sodium 
nitrite in normal volunteers during a 48-h intravenous infu- 
sion. 12 This study revealed that nitrite PK after short-term and 
long-term infusion was different. Nitrite systemic exposure in 
plasma increased less-than-proportionally at higher doses, 
which indicated a nonlinear PK. Methemoglobin concentra- 
tions also increased less-than-proportionally with increasing 
sodium nitrite doses, which suggested a nonlinear nitrite- 
exposure-toxicity relationship. The observed nonlinearity has 
not been well characterized and along with high variability 
observed in nitrite exposure and toxicity between subjects, 
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Figure 1 Nitrite-hemoglobin interactions in RBC. 

it poses a great challenge to appropriately select therapeutic 
doses for the patient population. 

To this end, we developed a multiscale systems pharma- 
cology model to broaden the understanding of such com- 
plexity and provide a quantitative pharmacological guide to 
determine optimal dosing regimens. Our model adequately 
characterizes NO metabolome and methemoglobin con- 
centrations at eight different sodium nitrite doses while 
considering the contributions of major components in the 
nitrite-hemoglobin system. The application of the developed 
model in sodium nitrite dosage selection may prove valuable 
in the clinical setting. 

RESULTS 

Normal volunteer data 

Concentrations of nitrite and nitrate in plasma and RBC, 
and methemoglobin were obtained from normal volunteers 
during and after a 48-h sodium nitrite infusion at escalating 
doses that ranged from 4.2 to 533.8 ug/kg/h (n = 12). 12 Data 
from one subject who received 533.8 ug/kg/h were excluded 
because of observed asymptomatic toxicity. Similarly, in 
another subject who received 445.7 ug/kg/h administration 
was stopped 3.2 h after infusion started because of toxicity 
and only available data were included. Nitrite and nitrate RBC 
concentrations of subjects who received doses from 4.2 to 
1 33.4 9 ug/kg/h (n = 6) were not available for analysis. A total 
of 333 plasma concentrations and 147 RBC concentrations 
of nitrite and nitrate, and 333 methemoglobin concentrations 
were used for model development. 

Structure model selection 

Typical representative individual concentration-time profiles 
of nitrite and nitrate in plasma and RBC as well as methemo- 
globin at maximal-tolerated dose (266.9 ug/h/kg) are shown 
in Figure 2. During the 48-h infusion, nitrite concentrations 
in plasma and RBC quickly reached a plateau within 6h, 
while nitrate slowly reached plateau by 12h. Methemoglobin 
rose slowly, reaching a maximum level within 12h. Only one 
subject who received 445.7 ug/kg/h showed asymptomatic 
methemoglobin concentrations above 5% immediately after 
the infusion ended. Plasma nitrite, plasma nitrate, and met- 
hemoglobin concentrations did not increase proportionally 
with increasing nitrite doses, which indicate dose-exposure- 
toxicity nonlinearity. Plasma and RBC nitrite concentrations 
rapidly returned to preinfusion concentrations. Plasma and 
RBC nitrate and methemoglobin concentrations remained 
elevated and increased immediately after the infusion ended 
before slowly returning to preinfusion concentrations. 



Individual time profiles of nitrite and nitrate in plasma and 
RBC and methemoglobin were fitted simultaneously to the 
proposed model (Figure 3). The developed model adequately 
described NO metabolome and methemoglobin disposition 
after a 48-h infusion of increasing nitrite doses. The final 
model parameters are summarized in Table 1. Key features 
of our final model include concentration-dependent increase 
in nitrite and nitrate clearance, nitrate entero-salivary recy- 
cling, and facilitated methemoglobin removal as described in 
detail below. 

Nitrite and nitrate disposition in plasma and tissue 

The less-than-proportional increase in nitrite and nitrate sys- 
temic exposure was explained by a concentration-dependent 
increase in clearance using a linear model. Increase in clear- 
ance was presumably due to saturable renal reabsorption at 
higher nitrite doses. 13 At maximal-tolerated dose, renal clear- 
ance was 1.159 l/min for nitrite and 0.014 l/min for nitrate, 
which was twofold and 1 .5-fold higher than renal clearances 
at baseline. Overall, the contribution of the renal excretion 
pathway of nitrite is small at endogenous levels as most of 
the excreted nitrite is reabsorbed (> 96%). 13 However, as 
renal reabsorption of nitrite is saturated at higher concen- 
trations, nitrite renal elimination becomes important and is 
needed in the final model to avoid nitrite systemic overes- 
timation. Nitrite metabolism by heme carriers, for example, 
myoglobin, CYP450, occurs in tissue under physiological 
conditions 14 and it was estimated to be 0.234 rr\\n-\k MYO ). 

Although nitrate is a nitrite metabolite subject to renal 
excretion, nitrate can be reduced back into nitrite via entero- 
salivary recirculation or reductase endogenously present in 
tissue. After dietary consumption (approximately 1-2 mmols 
per day), 15 nitrate is rapidly absorbed in the upper gastroin- 
testinal tract, accumulated in saliva, and reduced into nitrite 
by nitrate reductases in commensal microflora (e.g., Vionella 
spp). Some of the newly formed nitrite in saliva reaches the 
stomach and is reduced into NO or absorbed into the blood to 
interact with hemoglobin. 115 In our study, three subjects had 
high nitrate concentrations at predose (30.33 ±10.20 umol/l) 
and high plasma nitrite concentrations early after the initia- 
tion of infusion (Figure 4). Those subjects had higher nitrate 
value compared with the normal baseline (10. 78 ±1.33 
umol/l) to accommodate the dietary effect. We attempted to 
incorporate entero-salivary compartments that connected 
from plasma nitrate to plasma nitrite via additional transit 
compartments to represent physiological pathways in the gut 
and saliva. However, a simpler approach representing direct 
conversion of nitrate via a first-order process (k N02 RDT ) pro- 
vided similar improvement. Our model suggests that nitrate 
is reduced back into nitrite more readily via entero-salivary 
recirculation (1.65x10 -3 min -1 ) than in tissue by means of 
xanthine oxidoreductase (6.83 x10~ 5 min -1 ). 16 

Nitrite RBC uptake and nitrite-hemoglobin interaction 

Nitrite is rapidly sequestered by RBC and interacts to a greater 
extent with deoxyhemoglobin than oxyhemoglobin (1 .42x10~ 4 
vs.5.91 x10- 6 min- 1 -umoM), which contributes more extensively 
to methemoglobin formation, similar to deoxyhemoglobin/ 
oxyhemoglobin-nitrite interactions in vitro? 7 In addition, a 
small fraction (4.61%) of hemoglobin-nitrite-derived NO 
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Figure 2 Observed and model-predicted concentration-time profiles for nitrite and nitrate in plasma and RBC and methemoglobin after 
administration of maximal-tolerated dose (266.9 |jg/h/kg) for 48 h. 
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Figure 3 The quantitative systems pharmacology model of nitric oxide metabolome and methemoglobin. Dotted arrows and open boxes 
indicate stimulation process. Shaded compartments represent variables that were experimentally measured and open compartments 
represent variables with no measurements. 



contributes to methemoglobin formation via oxyhemoglobin, 
similar to observations in vitro that characterized NO con- 
sumption in RBC. 18 Our results suggest these interactions 
occur in vivo, which may contribute to greater NO production 
and NO-dependent vasodilation. 19 While nitrite-hemoglobin 
interactions are described with a second-order rate constant 



given hemoglobin abundance compared to nitrite, nitrite- 
deoxyhemoglobin (k N0 R -[deoxyHb] = 1.65min- 1 ), and nitrite- 
oxyhemoglobin interactions (k N03 R -[oxyHb] = 0.234 min -1 ) 
may behave as a pseudo-first-order process depending on 
nitrite's availability. However, a second-order rate constant 
can account for interindividual variability (I IV) in hemoglobin 
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Table 1 Summary of model parameter estimates. 


Parameter (Unit) 


Definition 


Population mean 


IIV (%CV) 


k PT N02 ( m ' n_1 ) 


Nitrite distribution rate constant for central-to-tissue CMT 


0.108 


9 


k TP N02 ( m ' n_1 ) 


Nitrite distribution rate constant for tissue-to-central CMT 


1.745 


140 


k PT N03 ( m ' n_1 ) 


Nitrate distribution rate constant for central-to-tissue CMT 


0.160 


- 


k TP N03 ( m ' n_1 ) 


Nitrate distribution rate constant for tissue-to-central CMT 


0.515 


- 


k N02 RDT ( m ' n_1 ) 


Nitrate-to-nitrite formation rate constant in central CMT 


1 .65E-03 


180 


k MYO ( m ' n_1 ) 


Nitrite-to-nitrate formation rate constant in tissue 


0.234 


- 


k T RDT ( m ' n_1 ) 


Nitrate-to-nitrite formation rate constant in tissue 


6.83E-05 


- 


CL (o L no2 C-m' 11 " 1 ) 


Basal clearance for endogenous nitrite 


0.382 


- 


C 

°N02 


Slope factor to nitrite clearance 


4.524 


- 


^N02 P 0 


Volume of distribution of nitrite in central CMT 


12.418 


82 


CL (0LNO3 ('- min " 1 ) 


Basal clearance for endogenous nitrate 


9.41 E-03 


150 


C 

°N03 


Slope factor to nitrate clearance 


0.017 


- 


V N0 3_P(') 


Volume of distribution of nitrate in central CMT 


12.418 


- 


k PR N02 ( m ' n_1 ) 


Nitrite uptake rate constant from central CMT into RBC 


0.019 


110 


k RP N03 ( m ' n_1 ) 


Nitrate diffusion rate constant from RBC into central CMT 


5.17E-03 


96 


Wr (min- 1 .umo|- 1 ) 


Formation rate constant of nitrate and methemoglobin from nitrite 


5.91 E-06 


46 


k No_R (min- 1 -umo|- 1 ) 


Formation rate constant of nitric oxide and methemoglobin from nitrite 


1 .42E-04 


- 


k DEG ( m ' n_1 ) 


Degradation rate constant for methemoglobin reductase 


0.016 


- 


q 

°MetHb 


Slope factor to methemoglobin reductase 


0.045 


140 


v R _ N02 (D 


Volume of distribution of nitrite in RBC 


5.816 


37 


v R _ N03 (i) 


Volume of distribution of nitrate in RBC 


4.391 


- 


^FLMethHb 0 


Volume of distribution of methemoglobin in RBC 


3.284 


50 


K HbNoi (min- 1 -umo|- 1 ) 


Formation rate constant of HbNO from nitric oxide in RBC 


2.86E-05 


- 


K HbN02 (min- 1 -umo|- 1 ) 


Formation rate constant of nitrate and methemoglobin from nitric oxide in RBC 


4.07-07 


- 


Residual variability 








Additive error for plasma nitrite (umol/l) 




0.024 




Proportional error for plasma nitrite (%) 




56.53 




Additive error for RBC nitrite (umol/l) 




0.217 




Proportional error for RBC nitrite (%) 




48.76 




Additive error for plasma and RBC nitrate (umol/l) 


3.848 




Proportional error for plasma and RBC nitrate (%) 


36.33 




Additive error for methemoglobin (umol/l) 


0.484 




Proportional error for methemoglobin (% 


) 


26.14 





CMT, Compartment. 



concentrations and allow characterizing ischemia in the 
patient population. 

Methemoglobin was formed less-than-proportionally at 
higher nitrite doses, which suggests either less formation or 
facilitated removal. 9 In our model, methemoglobin reductase 
activity (k MR(t) ) was induced to prevent methemoglobin accu- 
mulation. 20 The turnover of methemoglobin reductase activity 
was described by the indirect response model (IDR), 2122 and 
the enzyme activity was stimulated by methemoglobin in a 
linear fashion. The sensitivity term of methemoglobin reduc- 
tase activity (S MetHb ) may explain the unexpected increase in 
methemoglobin concentrations in one subject during infu- 
sion of a low dose (4.2 ug/kg/h). 12 The estimated S MetHb value 
for this subject was 73% lower than the population mean 
(0.045) resulting in less active methemoglobin reductase and 



consequently, higher methemoglobin concentrations. Simi- 
larly, one subject receiving 445.7 ug/kg/h showed higher- 
than-expected methemoglobin concentrations and a S MetHb 
term 88% lower than the population mean. Alternatively, 
subjects with high S MetHb values have higher sensitivity to 
activate methemoglobin reductase and are likely to produce 
less methemoglobin as seen in subjects who received doses 
from 8.3 to 66.8 ug/kg/h and whose S MetHb was twice to four 
times higher than the population mean. Overall, S MetHb values 
ranged from 0.0054 to 0.171, which were reflected in the 
high IIV for S MetHb (> 100%). We also tested limiting methe- 
moglobin formation with an upper tolerable-bound serving 
as a physiological limit. 23 However, facilitating the removal of 
methemoglobin improved the overall fitting and was therefore 
incorporated into the final model. 
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Figure 4 Observed and model predicted concentration-time profiles 
for plasma nitrate (a) and nitrite (b) in a subject with evidence of the 
entero-salivary circulation. 

Covariate analysis 

Potential relationships between model parameters (i.e., 
nitrite distribution to tissue compartment, nitrate formation 
in plasma, clearance of nitrate and methemoglobin) and 
relevant covariates (body weight, gender, and age) on NO 
metabolome disposition were tested; however, no statistically 
significant covariate was identified. We observed relation- 
ships among covariates (a negative correlation between age 
and nitrite infusion rate, and higher infusion rate for females 
than for males) which appeared to be an unintended result 
of the study design and were carefully considered during 
covariate analysis. 

Model evaluation 

Goodness-of-fit plots (Figure 5) show that model predictions 
were in close agreement with the observed NO metabolome 
and methemoglobin concentrations. Population and individ- 
ual predicted concentrations showed adequate distribution 
around the unity line (Figure 5a, b). However, model-predicted 
population data were slightly overestimated for plasma nitrite 
and high RBC nitrate and methemoglobin concentrations 
observed immediately after infusion termination were under- 
estimated. The individual weighted residuals did not reflect 
systemic deviations (Figure 5c,d). The prediction-corrected 
visual predictive checks revealed that the final model predic- 
tions were in a reasonable agreement with the observed val- 
ues (Supplementary Figure S1 online). As the covariance 
step did not converge successfully, no standard errors and 



correlation matrix could be obtained, which underscores the 
degree of remaining uncertainty regarding the parameter 
estimates. The bootstrapping was not practically feasible due 
to computational limitations. We used the REPEAT option in 
NONMEM to assess the reliability of the final estimates and 
identified k pT N03 (79%), k TP N03 (80%), and k N02 RDT (72%) to 
be the most variable parameters. Future studies with more 
data may be necessary to corroborate our findings. 

DISCUSSION 

We developed a physiological, systems pharmacology- 
based model to describe the disposition of nitrite, nitrate, 
and methemoglobin in plasma and RBC. The model explains 
nitrite-hemoglobin interactions and nonlinearity in dose- 
exposure-toxicity relationships after a long-term sodium 
nitrite infusion to healthy volunteers. 

Previous studies have reported renal excretion as the main 
elimination pathway for NO metabolic species. 5 ' 1324 Renal 
clearance of nitrate/nitrite in dogs was dose dependent after 
sodium nitrate administration due to increases in nitrate/ 
nitrite urinary excretion as nitrate doses increased from 0 to 
40 ug-min~ 1 -kg~ 1 while the amount reabsorbed back into cir- 
culation progressively decreased from 97 to 84. 6%. 13 In pre- 
vious studies that assessed nitrite PK in normal volunteers, 
dose dependency was not observed since high doses were 
administered and different dosages were not evaluated sepa- 
rately. 89 We observed lower systemic exposure of nitrite than 
expected, which suggested accelerated elimination, similar 
to previous findings in normal volunteers with higher than 
expected nitrate urinary excretion after nitrate ingestion. 524 
The CL R N03 (0.010 l-min -1 ) 5 was in close agreement with the 
CL R N03 we estimated at similar doses (0.013 l-min -1 ). We 
attempted to use a physiologically-based model to describe 
nitrite and nitrate clearances and define renal clearance 

aSCL R_NO x = CL { 0 )N O x ^- FR NO^ Where C W ePreSentS 

nitrite/nitrate filtration and tubular secretion before reabsorp- 
tion. Then, fractional reabsorption (FR ) would decrease 
from its basal tubular reabsorption (97%), 13 and clearance 
would increase with higher nitrite doses. Although our study 
has a wide range of doses, we were unable to use this 
approach or the Michael is-Menten equation due to the lim- 
ited number of subjects, and instead used a simpler model. 

Hon et al. reported that incorporating nitrate-to-nitrite con- 
version pathway did not show overall improvement, suggest- 
ing this is a relatively insignificant metabolic pathway during 
short-term nitrite infusion. 9 However, a significant improve- 
ment in overall fitting and reduction in objective function value 
(OFV) (dOFV = -125.81) was observed by including this 
pathway after long-term infusion. The necessity of nitrate-to- 
nitrite conversion was also evident from three subjects whose 
initial plasma nitrite and nitrate concentrations were higher 
than expected and from experiments using antiseptic mouth- 
wash in rats 25 and humans. 26 In humans, rinsing the mouth 
with antibacterial mouthwash prior to nitrate load abolished 
plasma nitrite rise up to 80% without changing nitrate accu- 
mulation in saliva, which suggests that nitrate-to-nitrite recy- 
cling pathways play an important role as recurring plasma 
nitrite sources 27-29 during prolonged nitrite administration. 
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Figure 5 Goodness of fit plots for population and individual predictions including observed vs. population-predicted (a) and individual- 
predicted concentrations (b); and individual weighted residuals vs. individual predicted concentrations (c), or time (d). Solid lines represent 
the line of unity. 



This observation suggests that nitrate-rich food consumption 
during treatment with nitrite should be controlled. 

The previous study with two short-term nitrite infusions 
showed that methemoglobin was lower during the second 
infusion than during the first infusion, which was modeled by 
two methemoglobin formation rates (0.155 vs. 0.090 min -1 ). 9 
Methemoglobin formation is tightly regulated to maintain 
homeostasis, 23 and its excessive accumulation is prevented 
by increased methemoglobin reductase activity, 20 which 
our model incorporated to describe the less-than-propor- 
tional increase in methemoglobin with increasing nitrite 
doses. k MR(t) and S MetHb explained the higher than expected 
increases in methemoglobin in two subjects. Individuals with 
lower enzyme activity (k MR(t) ) or with less enzyme induction 
capability (S MetHb ) may experience higher methemoglobin 
increase than expected. Concentrations of methemoglobin 
reductase may be measured in future studies to corroborate 
these results. 

Nitrite-hemoglobin interaction mechanisms are rather 
complex and have not been well established in vivo. Various 
hypotheses have tried to describe nitrite uptake into RBC. For 
instance, nitrite uptake into RBC may be the result of carrier- 
mediated nitrite anion (NO~) uptake 11 or nitrous acid diffu- 
sion. 1030 These hypotheses were tested in our model but no 
improvements were observed in the accuracy of model pre- 
dictability. Similarly, nitrite-hemoglobin and NO-hemoglobin 



interactions result in a number of intermediate reactions 
that have been widely studied in vitro, 17,1 8,31,32 but addition of 
intermediate steps did not show improvements in the model, 
similar to previous findings in rats that represented full cycles 
as a single step. 33 Additionally, nitrite-myoglobin interactions 
need to be further investigated in vivo. 34 

Several RBC NO release mechanisms have been pro- 
posed. Nitrite compartmentalization near the RBC membrane 
may facilitate a localized reaction between deoxyhemoglobin 
and nitrite for NO release 2 with a diffusional barrier formed 
by methemoglobin. 11 A second model suggests the forma- 
tion and diffusion of a less-reactive intermediate nitrosating 
agent such as dinitrogen trioxide. 3536 Other mechanisms 
involve S-nitrosothiol or HbNO formation which preserve and 
transport NO 1837-39 and escape RBC easily. 40 It is possible 
that NO uses various mechanisms to escape RBC, reach 
vascular endothelium, and exert its vasodilatory effect. 41 The 
next important step would be to identify which NO species 
would correlate with the observed vasodilatory effect and to 
incorporate such mechanism to build a comprehensive PK/ 
PD model that could be used in optimal therapeutic dosing 
selection. 

A limitation of our study was the small number of subjects 
overall and per-dose, which prevented us from exploring 
potential models that might better explain the observed data. 
We were unable to incorporate intermediate reactions steps, 
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which might be responsible for the abrupt peaks in some 
nitrite and nitrate concentration-time profiles. The selec- 
tion of the final model was on the basis of finding a balance 
between mechanistic detail and the parsimonious nature of 
modeling, given the data. 42 This could partially explain the 
degree of remaining residual error and uncertainty of some 
parameter estimates, which underlines the shortcomings of 
the present analysis and the need for more work to be done 
in the future. 

In summary, this study represents a quantitative systems 
pharmacology model for NO metabolome characterizing the 
time profiles of nitrite, nitrate and methemoglobin following 
a long-term intravenous infusion of sodium nitrite in healthy 
volunteers. Our study suggests that nitrate-nitrite recircula- 
tion via the entero-salivary pathway may influence nitrite con- 
centrations and needs to be controlled. As methemoglobin 
serves as a toxicity measure, pre-evaluation of methemoglo- 
bin reductase activity may be useful to assess tolerability to 
methemoglobin rise upon nitrite administration. Inclusion of 
major nitrite-related hemoglobin species could be important 
to account for varying concentrations of hemoglobin oxygen- 
ation status, especially when implementing to patients with 
low oxygenation conditions (e.g., ischemia). The present 
model provides important insights into the nitrite-hemoglobin 
interactions and the nonlinear disposition of nitrite, nitrate 
and methemoglobin in humans. The application of our model 
in sodium nitrite dosage selection for patient population will 
be valuable in clinical development. 

METHODS 

Study design. The data were obtained from a Phase I dose 
escalation study. 12 The study design, informed consent and 
ethics procedures, and clinical results have been reported 
elsewhere. 12 Twelve healthy individuals aged 21-56 years old 
(39 ±9), weighing 49-1 15 kg (77.8 ±19) were intravenously 
infused with escalating sodium nitrite doses at 4.2 (n= 1), 8.3 
(n= 1), 16.7 (n= 1), 33.4 (n= 1), 66.8 (n= 1), 133.4 (n= 1), 
266.9 (n = 3, maximal-tolerated dose), 445.7 (n = 2), and 
533.8 (n = 1) ug/h/kg for 48 h (Supplementary Figure S2 
online). Concentrations of nitrite, nitrate and methemoglobin 
in plasma and RBC were measured using chemilumines- 
cence and gas chromatography-mass spectrometry methods 
previously described. 12 

Pharmacokinetic analysis. Population PK model develop- 
ment and simulations were performed using NONMEM ver. 
7.2. (ICON Development Solutions, Ellicott City, MD). NON- 
MEM was compiled by Intel Visual Fortram Compiler 11.0 
compiler (Intel Corporation, Santa Clara, CA) on a Window 
Server 2008 R2 operating system. NONMEM PREDPP 
library subroutines ADVAN8 with TOL 4 were used. Nitrite, 
nitrate, and methemoglobin in plasma and RBC were fitted 
simultaneously to the proposed model shown in Figure 3. 
The model included distribution and interconversion in tis- 
sue, uptake into RBC, and nitrite-hemoglobin interactions for 
methemoglobin and RBC nitrate formation. The model devel- 
opment process is presented in Supplementary Table S1 
online. 



Nitrite and nitrate disposition in plasma and tissue. Upon 
nitrite administration, nitrite in plasma diffuses into RBC (k PR 
N02 ), distributes into tissue (k PTN02 and k TP N02 ), or is directly 
eliminated (CL R N02 ). Nitrate in plasma distributes into tissue 
( k PT_N03 and k TP~No 3 ) or is directly eliminated (CL R N03 ). Nitrite 
in tissue may undergo oxidation (k MY0 ) by heme-carriers such 
as myoglobin. The differential equations for plasma and tis- 
sue nitrite and nitrate are: 

1 df Jp = k in NQ2 + k TP N02 [NO~] T + k N02 _ RDT [NO~] p 

- k PT N02 [NO~] p - k PR NQ2 [NO~] p ( 1 ) 

-CL RN02{t) /V N02 _ p [NO-] p 



1 df Jr = k PT _ N02 [NO-] p + k TRDT [NO-\ 
-k TP N02 [NO-] T -k MYO [NO-\ 



(2) 



d[NO-] p 
dt 



[NO-] R + [NO~] T 

N0 2 _RDT 

[NO~] P 

k PT_N0 3 

[NO-] p (3) 

-OL RN03{t) /V N03 _ p -[NO-] P 



= k PT N03 ■ [NO~] P + k MY0 ■ [NO~] T 
dt ~ (4) 

-k TP _ NO -[N0 3 ] T -k T _ RDT -[N0 3 ] T 



where [NO~ ] p and [NO~ ] p represent nitrite and nitrate 
amount in central compartment, and [NO~ ] T and [NO~] T 
represent nitrite and nitrate amount in tissue. V N02 p and 
V NQ3 p represent the volumes of distribution for nitrite and 
nitrate in the central compartment. 

The increase in nitrite and nitrate clearances with increas- 
ing concentrations 12 was described with a linear slope fac- 
tor, S N02 or S N03 , in which renal clearance (CL R _ N02(t) and 
CL R N03(t) ) increased from their basal clearance (CL (0) _ NO2 
and"CL (0) N03 ) at endogenous concentrations, depending 
on the change in nitrite (A[NO~] p ) or nitrate (A[NO~] p ) 
concentrations: 



CL 



■R_N02(t) 



CL,, 



0)_NO 2 



(1 + S N0 -A[A/0 2 -] P ) 



CL 



■R_N03(t) 



= CL {0) NO3 -0 + S NO ^A[NO-] p 



(5) 



(6) 



As nitrite is endogenously present, 43 k jn N02 represents a 
zero-order production rate of endogenous nitrite, defined as 
a secondary parameter: 



x in_N0 2 rx el_N0 2 



lNO- 2 ] Pi +k el _ NO -[NO;l 



[NO] Rm [deoxyHb] 



(7) 



(0) 
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where k ej N02 and k el N03 correspond to first-order elimina- 
tion rate constant of nitrite and nitrate at baseline and are 
defined as CL (0) _ N0 /V N02 ,_ p and CL (0)N JV NO3P 

In our data, three out of 12 subjects started with higher 
predose nitrate concentrations compared with baseline and 
higher plasma nitrite concentrations at the early phase of the 
infusion. As nitrate converts back to nitrite via entero-salivary 
recycling by means of nitrate reductase, 124 we incorporated a 
first-order rate constant (k N02 RDT ). Similarly, tissue nitrate can 
also be converted back into nitrite (k T RDT ). 43 

Nitrite and nitrate in RBC. Once nitrite in plasma diffuses into 
RBC (k PR N02 ), it interacts with deoxyhemoglobin (deoxyHb) 
to form NO and methemoglobin (MetHb), and with oxyhe- 
moglobin (oxyHb) to form nitrate and methemoglobin, which 
was depicted by the second-order rate constants k N0 R and 
k NQ3 R . Nitrite-derived NO interacts with deoxyhemoglobin 
and oxyhemoglobin to form HbNO (k HbN01 ) or methemoglobin 
and nitrate (k HbN02 ). The nitrate formed in RBC diffuses back 
into plasma (k RP ^.The differential equations for RBC nitrite 
and nitrate are: 



d[NQ-] R 



■ rx PR_N0 2 

[NO-] P -k NO R [NO-] R 



• [deoxyHb] - k NOz _ R ■ [NO;] R ■ [oxyHb] 



d\k 

JL ^ M = Ksvn ■ (1 + S Hb , ■ A[MetHb] I V R Met 

' ^DEG ' [^/Wfi](f) 



where V R MetHb is methemoglobin volume of distribu- 
tion in RBC. k SYN was a secondary parameter defined 

Oxyhemoglobin and deoxyhemoglobin concentrations 
were defined as a fraction of the total known hemoglobin 
at predose from which methemoglobin concentrations has 
been subtracted. In normoxic conditions, hemoglobin is found 
prevalently in the oxygenated form (77%) and deoxygenated 
hemoglobin corresponds to 23%. 45 It was assumed that total 
hemoglobin remains constant given the long RBC lifespan. 
Therefore, the changes of oxyhemoglobin and deoxyhemo- 
globin were defined as: 

[oxyHb] {t) = ([Hb] total - [MetHb] {t) )- 0.77 (1 3) 

[deoxyHb\ t) = ([Hb] total - [MetHb] {t) )- 0.23 (1 4) 



The initial conditions are defined by the steady-state 
values: 



d[NO-] R 
dt 



d[NO] R 



= k N o 3 _ R • [NO~] R • [oxyHb] - k RPNOs 
■[NO- 3 ] R + [NO] R [oxyHb] 

= k No_R ■ \. no I\r ■ [deoxyHb] - k HbN02 



0) 



^y^, ivu_h l ^ J/-f l j j no/vu^ (10) 

• [NO] R • [oxyHb] - k HbNm ■ [NO] R ■ [deoxyHb] 

where [NO~ ] R , [NO~ ] R , and [NO] R represent nitrite, 
nitrate, and NO amount in RBC. 

Methemoglobin formation and loss. Methemoglobin results 
from interactions between nitrite and oxyhemoglobin 
(k N03 R ) or deoxyhemoglobin (k N0 R ) and between NO and 
oxyhemoglobin (k HbNQ2 ) and is converted back to oxygen- 
binding hemoglobin forms by means of methemoglobin 
reductase (k MR(t j). 44 We assumed that methemoglobin 
reductase activity (k MR ) is induced by methemoglobin 
rise 20 to characterize the less-than-proportional increase 
in methemoglobin at higher nitrite doses, since lack of 
enzyme activity is known to result in disproportionate met- 
hemoglobin accumulation. 44 The enzyme activity turnover 
was described with its production (k SYN ) and loss (/c D£G ), 21 ' 22 
and the enzyme production was assumed to be stimulated 
by the change in methemoglobin amount (A[MetHb]) with 
the slope factor S MetHb . 



d[MetHb] 
dt 



k No_R • [NO;] R • [deoxyHb] + k NOs R 

■ [NO-] R ■ [oxyHb] - [k MR \ t) ■ [MetHb] (11) 
+ k HbN02 -[NO] R -[oxyHb] 



k N0 _ R ■ [deoxyHb\ 0) + k NO R • [oxyHb] 



(15a) 



(0) 



k p T no [NO:] p +k TRDT -[N0 3 ] T 
[NO'], = PT - N ° 2 \ 2JP(0) t TRDT (1 5b) 



k TP_N0 2 + ^MYO 



[NO; 



k N o 3 _p ' [N0 2 ]p m + k TP WOa ■ [NO~ 3 ] Tm 
+k RP _ NO -[NO- 3 ]„ (15c) 
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[NO; 



D 



k . k 
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[NO; 
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[N0] = k NOR -[NO-] Rm -[deoxyHb] m ^ 

[deoxyHb\ 0) + 

^HbNO^ 

[oxyHb\ 0) 

The rate of methemoglobin reductase activity at baseline 
was a secondary parameter defined as: 

k N0R ■ [NO-] R{q) ■ [deoxyHb] {0) + k NOz _ R ■ [NO;]^ 

= [^y^] ( o) + k H bN0 • fr ■ [NO] R{0) • [oxyHb] {0) 
MRw ~ [MetHb\ 0) 

Methemoglobin [MetHb] {0) and [NO~] p baselines were 
obtained from the individual average of known predose mea- 
surements and were read directly from the data file which is 
available in Supplementary Data online. 

The final model included 24 primary parameters (Table 1) 
and three secondary parameters (k SYN , k MR(0) , and k jn NQ2 ) for 
five dependent variables (plasma nitrite and nitrate, RBC 
nitrite and nitrate, and methemoglobin) fitted simultaneously. 
Five parameters (k N02RDV k TRDV k N03 R , k HbN01 , and k NO R ) 
were redefined using a dummy variable to bring parameter 
estimation value close to 1. No covariance step was per- 
formed in the model. 

Random variability models.The exponential model was used for 
interindividual (I IV) variability. Initially, I IV terms were estimated 
for all model parameters and then removed from parameters 
with high r|-shrinkage (> 60%). The residual error for plasma 
nitrite and nitrate, RBC nitrite and nitrate, and methemoglo- 
bin was evaluated using a combination of the proportional and 
additive models. The model parameter estimation used first- 
order conditional estimation method with INTERACTION. 

Covariate analysis. Relevant covariates (body weight, gen- 
der, and age) were screened via forward selection and back- 
ward elimination steps. Covariate significance was evaluated 
based on changes in OFV, reduction in I IV, and overall model 
performance. A difference between the OFV greater than 3.84 
(P< 0.05) was considered significant and covariates included 
in the corresponding parameter. The goodness-of-fit for mod- 
els were assessed numerically (OFV, AIC, ETA shrinkage and 
so on) and graphically. Graphical diagnostics using Xpose 4.0 
(Uppsala University, Uppsala, Sweden) included scatter plots 
of observed, population, and individual predicted (PRED and 
IPRED) concentrations vs. time; observed vs. PRED; weighted 
residuals vs. time or IPRED; and ETA matrix correlation. 

The stability of the model was evaluated by prediction 
corrected visual predictive checks 46 in which parameter 
estimates were used to simulate the data for 1 ,000 virtual 
patients. The 5th, 25th, 50th (median), 75th, and 95th per- 
centiles were calculated and the distribution of simulated 
concentrations was visually compared with the measured 
concentrations at each sampling time point. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE OF THE 
TOPIC? 

S Nitrite, as an NO donor, has been shown to 
be safe for human administration. Following 
prolonged infusion, a nonlinear dose-exposure- 
toxicity relationship was described, which poses 
a challenge in further clinical development. 

WHAT QUESTION THIS STUDY ADDRESSED? 

y A physiological and mechanistic model was 
developed to explain nitrite-hemoglobin inter- 
actions and nonlinear dose-exposure-toxicity 
relationships in NO metabolome and methemo- 
globin after nitrite administration. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

S Dose-dependent increase in nitrite/nitrate clear- 
ance and methemoglobin removal at higher 
nitrite doses explain the less-than-proportional 
increase in exposure and toxicity. Nitrate-to- 
nitrite reduction occurs in tissue and via entero- 
salivary recirculation and influences nitrite PK. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

S Nitrite-nitrate recirculation may be relevant for 
nitrate-rich food intake and may require control 
during nitrite administration. As methemoglobin 
serves as a toxicity measure, pre-evaluation of 
methemoglobin reductase activity may deter- 
mine individuals' tolerability to methemoglobin 
rise upon nitrite administration. Nitrite-related 
hemoglobin species could be important in low- 
oxygenation conditions. The application of our 
findings in sodium nitrite dosage selection for 
patient population may prove valuable. 
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